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Abstract 

We carry out a state-of-the-art assessment of long baseline neutrino oscillation experiments 
with wide band beams. We describe the feasibility of an experimental program using existing 
high energy accelerator facilities, a new intense wide band neutrino beam (0-6 GeV) and a 
proposed large detector in a deep underground laboratory. We find that a decade-long program 
with 1 MW operation in the neutrino mode and 2 MW operation in the antineutrino mode, a 
baseline as long as the distance between Fermilab and the Homestake mine (1300 km) or the 
Henderson mine (1500 km), and a water Cherenkov detector with fiducial mass ~ 300 kT has 
optimum sensitivity to a nonzero #13, the mass hierarchy and to neutrino CP violation at the 
3cr C. L. for sin 2 2#i3 > 0.008. This program is capable of breaking the eight-fold degeneracy 
down to the octant degeneracy without additional external input. 
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1 Introduction 



There is abundant evidence that neutrinos oscillate and hence have mass . Atmospheric neutrinos 
oscillate with a mass-squared difference lAm 2 .^ ~ 0.0025 eV 2 and mixing angle #23 ~ 7 r/4 [2] while 
the corresponding parameters for solar neutrinos are Am^ ~8x 10~ 5 eV 2 and 9 U ~ n/6 0. The 
unknowns are the angle #13 coupling the atmospheric and solar sectors, the CP-violating phase <5cp> 
and the sign of Am^, which fixes the hierarchy of neutrino masses. All we presently know about #13 
is that sin 2 #13 < 0.05 at the la C. L. [3]. Measurements of the unknown parameters are the main 
goals of future long-baseline neutrino experiments. Complete and precise experimental knowledge 
of neutrino parameters is essential to test neutrino mass models. Moreover, a determination of the 
type of mass hierarchy is of importance for the feasibility of observing neutrinoless double beta if 
the heaviest neutrino is O(0.05) eV. An observation of this decay will provide much confidence in 
the Majorana nature of neutrinos. 

In the near future, experiments with low energy conventional neutrino beams will attempt to 
detect — > u e transitions at the atmospheric scale. Concurrently, reactor experiments will try 
to determine or constrain #13 by observing D e disappearance at the atmospheric scale. However, 
the complementary information gathered from both types of experiments will be inadequate to 
determine the mass hierarchy or to establish CP-violation in the lepton sector. In recognition of 
this fact, planning is underway for future facilities that will address these outstanding issues. 

Two types of neutrino beams are being considered for future long-baseline experiments: off-axis 
beams with narrow band energy spectra (as will be used in the Tokai-to-Kamioka (T2K) experiment 
at Japan Proton Accelerator Proton Complex (J-PARC) jS], and proposed for the NuMI Off-axis v e 
Appearance Experiment or NO^A and a wide band beam [7J|S| for which the neutrino energy 
is inferred from quasielastic scattering events. 

Off-axis beams jH] are advocated because their beam energy is narrow and depend to first 
order only on the kinematics of pion decay, so that a counting experiment can be carried out 
at an oscillation maximum. Another positive feature of the narrow beam is the reduction of a 
high energy tail which could contribute background to reconstructed events at low energy via 
neutral-current feed down. The features that make off-axis beams attractive also lead to their 
limitations. The very nature of a counting experiment precludes the possibility of using spectral 
energy information. Consequently, it is virtually impossible to resolve the eight-fold parameter 
degeneracy ^U] (which arises from three independent two-fold degeneracies (<!>cp,#i3), sgn(Am 2 1 ) 
and (623, ir / 2 — #23) E3^D)' unless multiple measurements with separate detectors are made at 
different off-axis angles and/or baselines. The capabilities and optimization of experiments with 
off-axis beams have been studied extensively |12j . 

Wide band beams possess the advantages of a higher on-axis flux and a broad energy spectrum. 
In an idealized sense, one may think of them as enabling many off-axis experiments in a single 
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experiment. Thus, parameter degeneracies may be resolved in a single wide band beam experiment. 
The disadvantage of a wide band beam is the need to determine the neutrino energy with good 
resolution and eliminate background from high energy tails of the spectrum. Therefore, detectors 
with good energy resolution and neutral-current rejection are needed to use the broad spectrum 
fruitfully and to reject the feed down from the high-energy tail efficiently 

The possibility of working with large detectors may not be too far into the future |13l 1141 . 
Two concrete proposals are under consideration for a Deep Underground Science and Engineering 
Laboratory (DUSEL) ^B] in the U.S., one at the Homestake mine in South Dakota, the other at 
the Henderson mine in Colorado. DUSEL will be designed to accomodate Mton class multipurpose 
detectors. The physics program of these detectors will include nucleon decay supernova neutrino 
detection, as well as long-baseline neutrino oscillation physics with an accelerator beam. 

The prospects of a very long-baseline experiment consisting of a wide band super neutrino beam 
and a large water Cherenkov detector were analysed in Refs. jHlEl) where, in the spirit of a first 
analysis, correlations and degeneracies were only partially accounted for. Preliminary work had 
suggested that sufficient background suppression could be achieved so as to extract the quasielastic 
events with sufficient purity. Since then, further detailed assessments with a full simulation of a 
water Cherenkov detector have been performed ^H]. The main result of that work is that neutral 
current backgrounds can be sufficiently suppressed while maintaining good signal efficiency with 
new software techniques. The remaining signal events are dominated by quasielastics, but have 
significant contamination of other v e induced charged current events. Nevertheless, neutrino energy 
resolution in the range of ~ 10% can be maintained. This work will be reported elsewhere in detail. 

In this paper, we perform a thorough state-of-the-art analysis with a realistic treatment of 
systematic errors, correlations and degeneracies so as to define an optimum program using a wide 
band beam. We consider accelerator and detector requirements for a feasible wide band experiment 
in the next two sections. In Section 4 we discuss our analysis methodology and present results of 
the analysis in Section 5. We summarize in Section 6. 

2 Accelerator Requirements 

A number of recent studies have examined possible intensities and neutrino spectra from U.S. 
proton accelerators j!91 1201 I21j . Here we summarize the understanding of the intensity versus 
proton energy available in the U.S. currently, in the near future, and with upgrades. 

High energy multi-MW proton beams are under consideration at FNAL. An examination is 
underway to increase the total power from the 120 GeV Main Injector (MI) complex after the 
Tevatron program ends |22| . In this scheme protons from the 8 GeV booster, operating at 15 Hz, 
will be stored in the antiproton accumulator (which becomes available after the shutdown of the 
Tevatron program) while the MI completes its acceleration cycle. Combining the techniques of 
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momentum stacking using the antiproton accumulator and slip-stacking using the recycler will 
raise the total intensity in the MI to > 1 MW at 120 GeV. In the ideal case, the length of the 
acceleration cycle is proportional to the proton energy, making the average beam power proportional 
to the final proton energy. However, there are fixed time intervals in the beginning and the end 
of the acceleration cycle for stable operation. These become important at low energies and reduce 
the performance below the ideal. Current projections suggest that ~ 0.5 MW operation between 
40 - 60 GeV and > 1 MW operation at 120 GeV is possible. 

Ambitious plans at FNAL call for a 8 GeV super-conducting LIN AC that can provide 1.5 x 10 14 
H~ ions at 10 Hz corresponding to 2 MW of total beam power |19| . Some of the 8 GeV ions could 
be injected into the MI to provide high proton beam power, 1 to 2 MW, at energy between 40 and 
120 GeV; e.g., 40 GeV at ~ 2 Hz or 120 GeV at ~ 0.67 Hz. Such a plan allows for flexibility in the 
choice of proton energy for neutrino production. 

The BNL Alternating Gradient Synchrotron (AGS) operating at 28 GeV currently can provide 
about 1/6 MW of beam power. This corresponds to an intensity of about 7 x 10 13 protons in a 2.5 
microsecond pulse every 2 seconds. The AGS complex can be upgraded to provide a total proton 
beam power of 1 MW |21| . The main components of the accelerator upgrade at BNL are a new 
1.2 GeV Superconducting LINAC to provide protons to the existing AGS, and new magnet power 
supplies to increase the ramp rate of the AGS magnetic field from about 0.5 Hz to 2.5 Hz. For 1 
MW operation the protons from the accelerator will be delivered in pulses of 9 x 10 13 protons at 2.5 
Hz. It has been determined that 2 MW operation of the AGS is also possible by further upgrading 
the synchrotron to 5 Hz repetition rate and with further modifications to the LINAC and the RF 
systems. A new neutrino beam could be built with conventional horn focussed technology and a 
200 m long pion decay tunnel placed on the slope of a specially built hill. Such a beam could be 
aimed at a detector ~ 2500 km away. 

To observe multiple oscillation nodes in the z/„ disappearance channel, it is necessary to have a 
wide band beam with energies from to 6 GeV. Protons with energy above ~ 20 GeV are needed 
to provide such a flux, which is clearly possible at both FNAL and BNL. The spectrum that is used 
in this paper was obtained using 28 GeV protons and a 200 meter long decay tunnel. For details 
of this spectrum see Refs. |17l I23j . Recent calculations have shown that similar, but more intense 
spectra per unit beam power can be obtained using 40 or 60 GeV protons from FNAL, where the 
meson decay tunnel could be made longer than 200 meters |24| . 

In our analysis we assume that the spectrum from either the FNAL or the BNL beam will be 
the same with a total average beam power of 1 MW. This allows us to make a proper comparison 
of the physics issues regarding the baselines. We comment on the additional flexibility obtained 
by the ability to change the proton energy at FNAL as well as the longer decay tunnel. We also 
comment on the impact on the sensitivity with 0.5 MW of average power from FNAL. It should be 
noted that 0.5 MW operation from FNAL does not appear to require a major upgrade. 
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If a large detector El US] is located at Homestake (HS) or Henderson (HD), the beam 
from FNAL (BNL) will have to traverse a distance of 1290 km or 1495 km (2540 km or 2770 km), 
respectively. At FNAL the inclination will be about 5.8° (HS) or 6.7° (HD). The existing experience 
at FNAL, from building the NuMI beam, could be extended to build a new beam to HS or HD. 
At BNL the beam would have to be built at an inclination of about 11.4° (HS) or 12.4° (HD). 
Current design for such a beam requires the construction of a hill with a height of about 50 m |21j . 
Such a hill will have the proton target at the top of the hill and a 200 m long decay tunnel on the 
downslope. In either case, it is adequate to have a decay tunnel with length shorter (200 to 400 m) 
than the NuMI tunnel (750 m) to obtain the needed flux. Since the focus is on lower energies, the 
beam intensity is maintained if the decay tunnel can be made about 4 m in diameter. 

Another advantage of the wide tunnel is that the option of running with a narrow band beam 
using the off-axis technique could be preserved. A 4 m diameter tunnel could permit the rotation 
of the target and horn assembly so that a 1° off-axis beam could be sent to the far detector. 

3 Detector Requirements 

A water Cherenkov detector with 300 kT fiducial mass can be built in the same manner as the 
SuperKamiokande detector (with 20 inch photo-multipliers placed on the inside detector surface 
covering approximately 40% of the total area) |25| by simply scaling it to larger size or by building 
several detector modules ^3 E]- Such a detector placed underground at DUSEL could have 
a low energy threshold (< 10 MeV), good energy resolution (~ 10%) for single particles, good 
muon/electron separation (<, 1%), and time resolution (<^ few ns). As noted earlier, it is important 
to obtain good energy resolution when using a wide band beam. This can be achieved in a water 
Cherenkov detector by separating quasielastic scattering events with well identified leptons in the 
final state from the rest of the charged-current (CC) events. The fraction of quasielastics in the 
total CC rate with the spectrum used in this paper is about 23% for the neutrino beam and 39% 
for the antineutrino beam. 

Separation of quasielastic events from the CC and neutral-current (NC) background is being 
used in the K2K experiment |26| . For the program considered here an essential problem is to 
separate electron shower events from other NC events, especially events containing a single 7r° in 
the final state. The goal is to search for v e induced showering events in the 0.5 to 4 GeV range. 
Single 7T° particles with energies of 1, 2, 3 and 4 GeV decay to two photons with a minimum 
and most probable opening angle of 16, 8, 5, and 4 degrees, respectively. The probability of a 
decay with an opening angle of more than 20° for 1, 2, 3 and 4 GeV 7r°'s is 40%, 8.2%, 3.6%, 
and 2.0%, respectively. In a water Cherenkov detector the position where the 7r° photons convert 
cannot be measured with sufficient precision from the pattern of Cherenkov light which tends to 
be two overlapping showering rings. At low 7r° energies the opening angle is sufficiently large 
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compared to the Cherenkov angle (42°) that single 7r°'s can be separated quite effectively. At 
energies greater than 2 GeV, however, the small angular separation between the two photons 
makes such separation difficult. It is well known that resonant single pion production in neutrino 
reactions has a rapidly falling cross section as a function of momentum transfer, q 2 , up to the 
kinematically allowed value |2Zj. This characteristic alone suppresses the background by more 
than 2 orders of magnitude for tt° (or shower) energies above 2 GeV |28j . Therefore a modest tt° 
background suppression (by a factor of ~ 15 below 2 GeV and ~ 2 above 2 GeV) makes the 7r° 
background manageable level over the entire spectrum. Such background suppression has recently 
been demonstrated using complete simulation and reconstruction in Ref. |18j . 

It has been argued that a liquid Argon time projection chamber (LARTPC) could be built with 
total mass approaching 100 kT [T^]. A fine grained detector such as this has excellent resolution for 
separating tracks, making it possible to use a large fraction of the CC events (rather than only the 
quasielastic events) to determine the neutrino energy spectrum. A LARTPC also has much better 
particle identification capability. Consequently, a LARTPC with a total fiducial mass of ~ 100 kT 
is expected to have similar performance as a 300 kT water Cherenkov detector. 

We assume a detector performance based on Ref. ^5j. For the physics sensitivities we assume 
1 MW operation for 5 years (at 1.7 x 10 7 sec/yr) in the neutrino mode, 2 MW operation for 5 
years in the antineutrino mode and a detector fiducial mass of 300 kT. With the running times, 
the accelerator power level, and the detector mass fixed, we consider baselines in the range of 500 
to 3000 km. With 1 MW of a 28 GeV proton beam, a baseline of 1300 km, and a 300kT fiducial 
volume detector we calculate ~ 230000 muon charged current and ~ 77000 neutral current events 
in 5 years of running in the neutrino mode in the absence of oscillations. Under the same running 
conditions in the antineutrino mode (with the horn current reversed) we find a total of ~ 74000 
antimuon charged current and ~ 27000 neutral current events; approximately 20% of the event rate 
in the antineutrino beam will be due to wrong-sign neutrino interactions. For both neutrino and 
antineutrino running, approximately ~ 0.7% of the CC rate will be from electron CC events which 
form a background to the Va — ► v e search. As an example, the total number of electron neutrino 
appearance events and the expected background for sin 2 29\3 = 0.1 and a normal hierarchy as a 
function of baseline are shown in Fig. ^ 

The event rates for other beam configurations using 40 or 60 GeV proton beam have been 
calculated in Ref. |24j . For equal proton beam power, higher proton energy could result in as much 
as a ~ 50% increase in the total event rate with a concomitant increase in the background. A 
careful evaluation of the optimization of event rate versus background could reduce the running 
time by a significant factor. Preliminary calculations suggest that the physics sensitivity using 
other beam configurations will remain approximately the same as the calculations reported in this 
paper after accounting for the effects of spectral variations on the backgrounds and signal. 
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4 Analysis techniques 



We simulate the expected event counts for the experimental program described above and compute 
physics sensitivities using the GLoBES software package j^H]- The description of the detector 
performance follows from the results in Ref. ^S] . We used this detector performance to create large 
numbers of Monte Carlo events for the following set of event classes, which in turn, are used to 
create migration matrices suitable for GLoBES: 

• v e appearance signal 

• disappearance signal 

• v e beam events 

• NC background mainly consisting of single tt° production. 

• Other NC events. 

• CC events with additional undetected particles. 

We have a set of six matrices for neutrinos and antineutrinos each. The background of antineutrinos 
in the neutrino beam is very small since it is suppressed both in production and detection by the 
cross section ratio of neutrinos to antineutrinos. For antineutrino running the opposite is true. 
There the neutrino (wrong-sign) contamination in the beam is nonnegligible due to the enhancement 
by the cross section ratio. Therefore, we account for this background via 6 additional matrices that 
correspond to the ones above for treating the neutrinos in the antineutrino beam. Our calculations 
show that the background due to tau neutrino interactions is low and can be safely ignored because 
the spectrum we use has low energy. When combined with the high threshold (3.5 GeV) for tau 
production as well as the low v T CC cross section, the background is not significant. If higher 
energy protons are used to create a spectrum with a high energy tail, then this background must 
be included. We have carefully checked that the GLoBES calculation of event rates agrees with the 
direct result of the Monte Carlo for the case of no oscillations. 

Once event rates are computed, the next step is the calculation of a \ 2 function and the 
inclusion of systematic errors. We use a \ 2 function for Poissonian processes as given, e.g., in 
Ref. [3U] . Systematic errors are implemented using a pull approach as in Ref. [3^ ■ The systematics 
considered are a normalization error for each type of signal, i.e., v e appearance, v e appearance, v„ 
disappearance and disappearance of 1% each, uncorrelated between the four types of signal. For 
the sum of all backgrounds to each signal we assume a 10% uncertainty in the normalization, again 
uncorrelated between the four types of signal. 
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As input or true values for the oscillation parameters and the Gaussian la ranges (in anticipation 
of precision measurements from near future experiments) we use: 

012 = 0.55 ± 10% , Amlx = (8.0 ± 0.8) • 10~ 5 eV 2 , 
6*23 = vr/4± 5%, Am§! = (2.5 ±0.125) • 10~ 3 eV 2 . (1) 

We include a 5% error on the matter density. Our analysis includes the correlations between all 
parameters and properly accounts for possible degeneracies. 

All calculations assume a normal mass hierarchy as input, whereas the fit always extends to 
the case of an inverted hierarchy. It is known that there are no qualitative differences in the 
sensitivities if an inverted hierarchy is assumed as the input; see e.g., Ref. jSl]. In that case the 
matter enhancement moves from neutrinos to antineutrinos, but since we assume nearly symmetric 
neutrino and antineutrino running, the results are not affected significantly. 

To allow a concise presentation of the physics results we define performance indicators for the 
various measurements. For 9±3 we choose the discovery of nonzero sin 2 20i3. Here a nonzero value 
of sin 2 20i3 is chosen as the true value in the simulation and a fit with sin 2 20i3 = is performed. 
The reach for discovering CP violation is computed by choosing a value for Sqp as input and fitting 
it with the two CP conserving values of 5qp = 0, ir. All other parameters (including $13) are free 
and the cases of an inverted mass hierarchy and 023 — ► 7r /2 — 023 are taken into account. For the 
exclusion of an inverted mass hierarchy a point in parameter space with normal mass hierarchy is 
chosen as the true value and the solution with the smallest x 2 value with inverted hierarchy has 
to be determined (global minimum of the x 2 function). The same applies to the resolution of the 
octant of 023- Data is generated with 023 < 7r/4 and the global minimum of the x 2 function with 
023 > vr/4 has to be found. 

The x 2 values obtained by the above procedure are converted to confidence levels by using the 
X 2 distribution for 1 degree of freedom for each of the above cases, i.e., \J Ax 2 corresponds to the 
significance in Gaussian standard deviations. 

5 Results 

Our results are: 

• Figure [21 shows the potential for discovering a nonzero value of 0i3 . The sensitivity does not 
depend strongly on the value of 5qp. A finite value of 0i3 can be established at 3a for sin 2 20i3 
as low as 0.005, irrespective of the value of 5qp- 

In what follows, we use the idea of a CP fraction, the definition of which is provided in Fig. [3] 
via illustration. For example, from the third panel one can see that for 75% of the values of 
Sep, the experiment can detect nonzero #13 at the 3a C. L. if sin 2 20i3 > 0.004. The most 
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optimistic case arises for a CP fraction of since the true value of the CP phase is the one 
for which the sensitivity is maximum. 

From Fig. 01 it is evident that the baseline is not crucial in determining the ability to confirm 
if #i3 is nonzero. This independence stems from the fact that with increasing baseline L the 
signal event rates do not drop as L~ 2 but much more slowly due to enhancement by matter 
effects. However, the background drops faster than L~ 2 since the oscillation with the solar 
mass splitting reduces the beam intrinsic background. This can be seen from Fig. |I] As a 
result, the ratio of signal to the square root of the background is approximately constant. 

In Fig. [5] we address the potential for discovering a normal mass hierarchy with a baseline 
of 1300 km. A normal hierarchy can be confirmed independently of 5c p at > 3<r so long as 
sin 2 20i3 is larger than about 0.01. For the inverted hierarchy the result is approximately the 
same because of the nearly symmetric v and v running. 

Figure El confirms the expectation that a longer baseline improves the sensitivity to the mass 
hierarchy since the matter effects increase with baseline. Note that an experiment with 
baseline below 1000 km has poor sensitivity to the mass hierarchy. However, the sensitivity 
(approximately) plateaus for baselines above 1500 km. As we see below, the optimal baseline 
is determined primarily by the discovery reach for CP violation. 

CP violation in the neutrino sector will be discovered if <5cp = and n can be excluded. 
Figure [7| shows the sensitivity for this measurement at two baselines, 730 km (left-hand panel) 
and 1300 km (right-hand panel). At a baseline of 730 km the mass hierarchy cannot be resolved 
and hence the sensitivity for 5cp > (for the inverted hierarchy, this occurs for Sep < 0) 
is severely limited by so-called 7r-transit: for certain combinations of true #13 and 5qp it is 
possible to fit the data with the wrong mass hierarchy and 5cp = vr |31j . The existence of 
this solution inhibits the ability to discover CP violation for a range of 5cp in an experiment 
with baseline from FNAL to the Soudan mine. At a baseline of 1300 km the mass hierarchy 
is well resolved (see Fig. |SJ) and therefore no problems arise due to 7r-transit. Baselines below 
approximately 1000 km are not optimal for studying CP violation. However, also at longer 
baselines it can be seen that the sensitivity to CP violation is greater for 5qp < 0. This 
is mainly due to the somewhat lower statistical significance of the antineutrino data; in the 
antineutrino mode, the signal event rate is lower and the background event rate is higher. 

In Fig. |SJ we show results for a CP fraction equal to 0.75 instead of for unity, since it is 
impossible for an experiment to have sensitivity to CP violation if 5qp = or ir. Note that the 
sensitivity decreases with baseline above 1000 km, which is to be attributed to the competition 
of genuine CP violation and fake CP violation induced by matter effects. Nevertheless, the 
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dependence of the sensitivity on the solar and atmospheric oscillation parameters and on the 
length of the baseline beyond 1000 km is weak and could be partly compensated by optimizing 
the low energy part of the spectrum. The flexibility at FNAL to vary the proton energy while 
keeping the total beam power constant is important for this optimization. Experiments with 
beams from FNAL to the Homestake mine or the Henderson mine, which have baselines of 
1300 km and 1500 km, respectively, have good sensitivity to CP violation compared to the 
optimal case of 1000 km and also have very good sensitivity to #13 and the mass hierarchy. 
Thus, baselines around 1300 - 1500 km represent a global optimum. 

• We now consider the prospect of breaking the (823,^/2 — 823) degeneracy, which is only 
present if the true value of 823 ^ ir/4. We display the results for two possible true values of 
823 in Fig. |5J To emphasize how challenging it is to break this degeneracy, we have selected 
values of 823 that are far outside the la range in Eq. (^Q). The ability to exclude the wrong 
octant depends on the choice of the true values of #13 and Sep- This dependence is not 
pronounced for most parts of the parameter space and therefore we chose to consider only 
the most pessimistic case, i.e. only that combination of true #13 and Sep which gives the 
lowest x 2 difference between the true and wrong octant. Hence, the results shown in Fig. [3] 
are conservative for the representative true values of 823- Clearly this measurement favors 
baselines longer than 2000 km. As noted earlier, such long baselines are not preferred for 
the detection of CP violation. Distinguishing the octant is a very difficult measurement 
which cannot be guaranteed to succeed with high confidence. If the neutrino spectrum at low 
energies can be enhanced, perhaps in a separate run with lower energy protons or off-axis 
geometry, then the signal to background ratio for this measurement could be improved. 

• Finally, we assess how our results depend on exposure and on the uncertainty in the overall 
normalization of the background. We have not considered the systematic uncertainties in 
the shape of the background assuming that the background shape, which is sharply peaked 
at low energies, will be measured sufficiently well by a near detector, as in the MINOS 
experiment |33| . 

From Fig. El we see that reducing the exposure by a factor of 3, degrades the sensitivity to 
a nonzero sin 2 2813 and to determine the mass hierarchy by about a factor of 2 in sin 2 2813, 
while the sensitivity to CP violation worsens by a factor of about 2 — 5 depending on the value 
of <5cp; if ^cp h as an unfavorable value, the experiment needs the highest possible exposure to 
maintain its ability to detect CP violation. Nevertheless, operation with half the total beam 
power, which could be achieved by the accumulator based upgrades at FNAL, is shown to 
have excellent sensitivity. 

Similar statements can be made for the dependence on the systematic uncertainty from 
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Fig. 1111 Again, the sensitivity to CP violation puts the most stringent requirements on 
the experiment. The most striking conclusion is that sensitivity to the mass hierarchy is only 
weakly affected by the systematic uncertainty. 

6 Summary 

We have studied the measurement of CP violation in the neutrino sector using a powerful, but 
conventional neutrino beam, using a MW-class proton source located in the U.S. (either at FNAL 
or BNL). The beam would be delivered to a large detector with fiducial mass ~ 300kT over a 
distance > 1000km. Such a detector could be built at a new Deep Underground Science and 
Engineering Laboratory (DUSEL) which is in the process of being evaluated in the U.S. The two 
possible sites for DUSEL at the Homestake mine in South Dakota and the Henderson mine in 
Colorado are at distances of 1290 km and 1487 km from FNAL, respectively (or 2540 km and 2770 
km from BNL, respectively). We have calculated the scientific reach and its dependence on total 
exposure and systematics for a range of distances that encompass these possibilities. 

The experiment is motivated by the need to have sensitivity to both the atmospheric and solar 
oscillation scales (a necessary condition for observing CP violation in the 3-generation framework), 
and to obtain an oscillatory pattern in the energy spectrum of neutrinos. Access to multiple nodes 
of oscillations leads to improved sensitivity for CP violation because the CP effect grows larger for 
successive nodes and because parameter ambiguities can be resolved. 

We have shown that the sensitivity to nonzero #13 using the appearance of electron neutrinos 
and antineutrinos is roughly independent of baseline in the range of 1000 to 3000 km. Also, 
the sensitivity to the mass hierarchy plateaus for baselines above 1500 km. On the other hand, 
sensitivity to CP violation decreases with baseline above 1000 km, but it is a weak function of the 
baseline beyond 1000 km. The baselines for a beam from FNAL to Homestake or Henderson are 
optimal for establishing a nonzero #13, the mass hierarchy and CP violation. The determination of 
the 623 octant, however, favors baselines longer than 1500 km. Both the 823 measurement and the 
CP measurement could be improved by increasing the low energy part of the flux. The sensitivities 
for a 1300 km baseline are shown in Figs. [21 El an d For the same baseline, the dependence 
of the sensitivities on exposure and systematic uncertainties are summarized in Figs. EH an d 111! 
respectively. 
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Figure 1: Event rates for neutrinos (left-hand panel) and antineutrinos (right-hand panel) as a 
function of baseline for sin 2 29\^ = 0.1 and a normal hierarchy. The bold lines show the signal for 
various choices of <5cp- The thin line shows the total background, which includes background from 
beam contamination and neutral current events. 




Figure 2: Discovery potential for sin 2 29i% ^ at a baseline of 1300 km. The bold iso-x 2 lines are 
3,4,5cr (from left to right) and the light lines show an increase of x 2 by 1. For all points to the 
right of the rightmost bold line, a nonzero value of sin 2 2#i3 can be established with at least 5a 
significance. 
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Figure 3: Schematic definition of the CP fraction. First, for each value of sin 2 2#i3, the result in 
the #13 — <5cp plane is used to determine the fraction of all CP phases for which there is sensitivity 
at the given confidence level (black arrows). Here we have taken the 3a curve of Fig. [3 Repeating 
this for all possible values of sin 2 26i% yields the figure shown in the middle panel. In the last step 
one chooses a set of values for the CP fraction and translates them into values of sin 2 26\^ by using 
the result from the second step. 
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Figure 4: Discovery reach for sin 2 2#i3 ^ at 3d for CP fractions (lowermost line, best case), 0.5 
(middle line) and 1 (uppermost line, worst case) as a function of the baseline. The detector mass, 
beam power and exposure are kept the same for all baselines. 
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Figure 5: Discovery potential for a normal mass hierarchy at a baseline of 1300 km. The bold 
iso-x 2 lines are 3,4, 5 a (from left to right) and the light lines show an increase of x 2 by 1. For the 
inverted hierarchy the results are approximately the same because of the approximately symmetric 
v and v running. 
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Figure 6: Discovery reach for a normal mass hierarchy at 3a for CP fractions (lowermost line, 
best case), 0.5 (middle line) and 1 (uppermost line, worst case) as a function of the baseline. The 
detector mass, beam power and exposure are kept the same for all baselines. 
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Figure 7: Discovery potential for CP violation at baselines of 730 km (left-hand panel) and 1300 km 
(right-hand panel). The bold iso-x 2 lines are 3,4, 5 a (from left to right) and the light lines show 
an increase of \ 2 by 1. For all points to the right of the rightmost bold line, CP violation can be 
established with at least 5cr significance. 
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Figure 8: Discovery reach for CP violation at 3<r for CP fractions (lowermost line, best case), 
0.5 (middle line) and 0.75 (uppermost line) as a function of the baseline. The detector mass, beam 
power and exposure are kept the same for all baselines. 
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Figure 9: Discovery reach for the octant of #23- Only the most conservative case with respect to 
the true values of 6*13 and <5cp is considered. The \ 2 difference between the true and wrong octant 
is shown as a function of the baseline for two representative true values of #23 that are far outside 
the la range in Eq. (^Q) (so as to emphasize how challenging this measurement is). The detector 
mass, beam power and exposure are kept the same for all baselines. 
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Figure 10: Dependence of the 3cr sensitivities on exposure for an experiment with 1300 km baseline. 
The two cases considered are 1/2 and 1/3 of the (full) exposure we have used throughout. 
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Figure 11: Dependence of the 3<r sensitivities on the uncertainty in the overall normalization of 
the background for an experiment with a 1300 km baseline. We have adopted a 10% uncertainty 
throughout. 
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